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Quasi-spherical subsonic accretion can be realized in slowly rotating wind-fed X-ray pulsars at 
moderate and small X-ray luminosities L x <4x 10 36 erg/s. In this regime the accreting matter set- 
tles down subsonically onto the rotating magnetosphere, forming an extended quasi-static shell. 
The shell mediates the angular momentum removal from the rotating neutron star magnetosphere 
by shear turbulent viscosity in the boundary layer near the magnetosphere or via large-scale con- 
vective motions. In the last case the differential rotation law in the shell is close to iso-angular- 
momentum rotation with (0 ~ 1 /R 2 . The accretion rate through the shell is determined by the 
ability of the plasma to enter the magnetosphere due to Rayleigh-Taylor instabilities while taking 
cooling into account. Measurements of spin-up/spin-down rates of quasi-spherically wind accret- 
ing X-ray pulsars in equilibrium with known orbital periods (like e.g. GX 301-2 and Vela X-l) 
enable determination the main dimensionless parameters of the model, as well as the estimate of 
the neutron star magnetic field. For equilibrium pulsars with independent measurements of the 
magnetic field, the velocity of the stellar wind from the companion can be estimated without the 
use of complicated spectroscopical measurements. For non-equilibrium pulsars, there exists a 
maximum possible spin-down torque exerted on the accreting neutron star. From observations of 
the spin-down rate and X-ray luminosity in such pulsars (GX 1+4, SXP 1062, 4U 2206+54, etc.) 
a lower limit on the neutron star magnetic field is derived, which in all cases turns out to be close 
to the standard one and in agreement with cyclotron line measurements. 
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1. Two regimes of quasi-spherical accretion in X-ray pulsars 

There can be two different regimes of quasi-spherical accretion (see e.g. §\§ for a recent review 
of previous studies of wind accretion). The captured stellar wind heated up in the bow shock at 
Bondi radius ~ Rb = 2GM/v 2 (where v is the relative stellar wind velocity) to high temperatures 
ksT ~ nipV 2 . If the characteristic cooling time of plasma t coo i is smaller than the free-fall time 
tff = Rb/ \/2GM/Rg, the gas falls supersonically toward the magnetosphere with the formation of 
a shock. This regime is usually considered in connection with bright XPSRs The role of X- 

ray photons generated near the NS surface is two-fold: first, they heat up plasma in the bow-shock 
zone via photoionization, and second, they cool down the hot plasma near the magnetosphere (with 
ksT ~ GM/Ra) by Compton processes thus allowing matter to enter the magnetosphere via the 
Rayleigh-Taylor instability 

In the free-fall accretion regime, the X-ray luminosity (the mass accretion rate M) is deter- 
mined by the rate of gravitational capture of stellar wind at Rb (Bondi-Hoyle-Littleton formula, 
M ~ pvR 2 B ). The accretion torque exerted on the NS due to plasma-magnetopshere interaction is 
always of the same sign as the specific angular momentum of the gravitationally captured stellar 
wind j m , and the NS can spin-up or spin down. 

If the relative wind velocity v at Rb is slow (< 80 km/s), the photoionization heating of plasma 
is important, but the radiation cooling time of plasma is shorter than the free-fall time, so the 
free-fall accretion regime is realized. If the wind velocity is larger than ~ 80 — 100 km/s, the 
post-shock temperature is higher than 5 x 10 5 K (the maximum photoionization temperature for 
a photon temperature of several keV); for L x < 4 x 10 36 erg/s, the plasma radiative cooling time 
is longer than the free-fall time, so a hot quasi-spherical shell is formed above the magnetosphere 
with temperature determined by hydrostatic equilibrium [||], [0]. Accretion of matter through such 
a shell is subsonic, so no shock is formed above the magnetosphere. The accretion rate M is now 
determined by the ability of hot plasma to enter magnetosphere. This is the settling accretion 
regime. 

2. Settling accretion regime: theory 

Theory of settling accretion regime was elaborated in [^J. In this regime, the accreting mat- 
ter subsonically settles down onto the rotating magnetosphere forming an extended quasi-static 
shell. This shell mediates the angular momentum transfer to/from the rotating NS magnetosphere 
by viscous stresses due to large-scale convective motions and turbulence. The settling regime of 
accretion can be realized for moderate accretion rates M < M* ~ 4 x 10 16 g/s. At higher accretion 
rates, a free-fall gap above the neutron star magnetosphere appears due to rapid Compton cooling, 
and accretion becomes highly non-stationary. 

Mass accretion rate through the hot shell is determined by mean velocity of matter entering 
the magnetosphere, u{Ra) = f(u)^2GM /Ra- The dimensionless factor f[u) is determined by the 
Compton cooling of plasma above magnetosphere and the critical temperature for Rayleigh-Taylor 
instability to develop and is found to be f(u) ~ (tff /t coo i) 1 / 3 . In the case of Compton cooling 

f(u)n0.4l*}i u n£ /n , (2.1) 
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where = M/[10 16 g/s] and = jx /[10 30 Gcm 3 ] is the NS magnetic moment. The definition of 
the Alfven radius in this case is different from the value used for disk accretion: 

/? A «1.6xl0 9 [cm](|^ 1 . (2.2) 

Accretion torques applied to NS in this regime are determined not only by the specific angular 
momentum of captured matter j m ~ 0.},R B (&b is the orbital angular velocity of the NS), as is the 
case of the free-fall accretion, but also by the possibility to transfer angular momentum to/from 
the rotating magnetosphere through the shell by large-scale convective motions. The plasma- 
magnetosphere interaction results in emerging of the toroidal magnetic field B t /B p = {K\ /Q (0^ — 
C0*)/(Ok{Ra), where co m is the angular frequency of matter at the Alfven radus, K\ ~ 1 the dimen- 
sionless coupling coefficient which is different in different sources, (Ok(Ra) is the Keplerian angular 
frequecy, and £ is the size of the region of plasma-magnetopshere angular moentum coupling in 
units of the Alfven radius R A . The NS spin evolution equation reads: 

Ico* = ^K 2 ^^^- +z MRl, (2.3) 
C R\ (Ok{Ra) 

where the second term takes into acount the angular momentum braught to the NS with the infalling 
matter (z < 1). This formula can be rearraged to the form 

la* = ZMR\{(o m - <o*) +zMR 2 A co* , (2.4) 

where the coupling coefficient is Z 036{K\ / QM^^ 11 pL^ 11 . 

The gas-dynamical treatment of the problem of angular momentum transfer through the shell 
by viscous turbulence stresses [||] showed that co m « ^(Ra/Rb)' 1 , where the index n depends on 
the character of turbulence in the shell. For example, in the case of isotropic near-sonic turbulence 
we obtain n ~ 3/2, i.e. quasi-Keplerian rotation distribution. In the more likely case of strongly 
anisotropic turbulence (because of strong convection) we find n^2, i.e. an iso-angular-momentum 
distribution 1 . 



3. Settling accretion regime: observations 

Equilibrium X-ray pulsars. In equilibrium XPSRs <&>*>= (e.g. Vela X-l and GX 
301-2). In this case, measurements of spin-up/spin-down near the equilibrium pulsar period P* q 
(or frequency co* q ) provides additional quantity dco* /dM (or dco*/dy, where y = MM eq is mass 
accretion rate or X-ray luminosity noprmalized to the equilibrium value). In this case (see 
[^] for more details) for the convective shell (n = 2) we find: a) equilibrium NS period via binary 
orbital period Pj,, mass accretion rate M (or X-ray luminosity L x = 0. IMc 2 ), NS magnetic field /I30 
and relative stellar wind velocity vg = v/ (lOOOkm/s) 

P* eq « \WO[s}pLll ll \P b /lO&)M- i *l n vt-, (3.1) 

1 If there is no convection in the shell, the magnetosphere interacts with the shell in a turbulent boundary layer. In 
that case the spin-down torque is ~ fi 2 /R^ (Rc = (GM/ft)* 2 ) 1 / 3 is the corotation radius) and is independent on M. This 
case of weak coupling can be realized for very faint slowly rotating XPSRs. 
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b) estimate of the coupling parameters Z eq or (K\/Q via P* and doo* /dy: 



j dm* I 
dM |e <? 



-eg 



— m*R 2 

11 w n A 



1. 



^fi)*/^y|y=l 



10- 12 rad/s 

c) estimate of the NS magnetic field via P* and dco* /dy 

'd(b*/dy\y=i 



(p*noos)M; 6 vn ^ n/n 



(3.2) 



A*30, 



eq 



10- 12 rad/s 
d) estimate of the stellar wind velocity 



(P*/100j) 



M 



-3/11 



16 



(3.3) 



; „,3 (1 - Z/Z «,)-M"^'(^) 



1/4 



(3.4) 



(note here the weak dependence on M and }x). The observed and calculated parameters of the 
equilibrium wind-fed pulsars Vela X-l and GX 301-2 are listed in Table 1. Note close values of 
the coupling parameter Z eq ~ 3 (or £ ~ 1/10) for both pulsars, and the independent measurement 
of the stellar wind velocity similar to the observed values. 

Non-equilibrium X-ray pulsars. From Eq. ( p3| ) it can be shown that CO* as a function of M 
reaches a mimimum at some accretion rate y cr = M cr /M eq = (3/2n + 3) 11 / 2 " < 1. For n = 2 we 
find: 



co, 



sd,max 



■1.13 x 10-H,sl(l- z /Z>V< (Z_)-" 4 (A)" 4 . (3 . 5) 



Aty < y cr the spin-down torque should anti-correlate with the the X-ray flux variations, deb* /dy < 
0, with ob s d ~ — R^ 3 ~ — M 6 l n . This is the case observed in long-term spinning-down XPSR GX 
1+4 [0, [|TTT]. From the condition \ob* d \ < \ob* d max \ we obtain the lower limit of the NS magnetic 
field: 



M30 > ^30 



:0.94 



CO' 



sd 



10- 12 rad/s 



-1/2 



-3/2 



100s 



7/8 



10d 



-3/8 



(3.6) 



If spin-up torque can be neglected, we find another estimate of the lower limit to the NS magnetic 
field 



M30 > /%) 



1.66 





11/13 ^ 


10- 12 rad/s 





'Ki 

c 



-11/13 



-4- M 



>- 3 /13 



16 



VIOOs 



11/13 



(3.7) 



Note that in contrast to Eq. (pTq), this estimate is independent of the poorly known stellar wind 
velocity vs and binary orbital period 

With decreasing M in non-equilibrium pulsars, the ratio of the toroidal to poloidal magnetic 



11/6 

M30 



2/3 



Below this 



field components increases, reaching B t ~ B p at M* 6 « 0.27 10 -i2 ra( j/ s 
luminosity acrretion becomes more non-stationary (likely the case of GX 1+4), but it is not still 
centrifugally prohibited. The propeller stage begins once Ra > R c at much smaller luminosities: 
M\l wO.OOSjuloCPVlOOs)- 11 / 3 . 

Very slowly rotating XPSRs There are known several very slowly rotating XPSRs, including 
some in HMXB (SXP 1062 with P* = 1062 s @, 4U 2206+54 with P* = 5550 s 0) and some 
in LMXB (e.g. 3A 1954+319, P* = 19400 s [fill]). Assuming disk accretion in such pulsars would 
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Table 1: Parameters of pulsars discussed 



Pulsars 



Equilibrium 



Non-equilibrium 



GX301-2 VelaX-1 



Ref. 



GX1+4 



SXP1062 4U2206 + 54 



11 



13 



Measured parameters 



P*(s) 


680 


283 


140 


1062 


5560 


ft(d) 


41.5 


8.96 


1161 




19 


v w (km/s) 


300 


700 


200 


~300* 


350 




2.7 


1.2 


? 


? 


1.7 


M 16 


3 


3 


1 


0.6 


0.2 


f | y =i(rad/s) 


1.5-10 


12 1.2- 10 12 


n/a 


n/a 


n/a 










-2.34 -10- 11 


-1.63 - 10 11 


-9.4- 10~ 14 


Obtained parameters 


/(") 


0.53 


0.57 








T 


14 


10 






> 8 


Z 


3.7 


2.6 








Bt/B p 


0.17 


0.22 








R A {cra) 


2 • 10' 


' 1.4- 10 9 








C0*/C0 K (R A ) 


0.07 


0.08 








fw,«»(W s ) 


500 


740 
















AC„«20 





' Estimated from the Corbet diagram 
* Typical velocity assumed in Be X-ray binaries 



require incredibly high NS magnetic fields (see e.g. discussion in [|15|]). However, application 
of our model to these and other non-equilibrium XPSRs (see Table 1) gives the low limits of the 
NS surface magnetic field in the usual range 10 12 — 10 G, and it is too preliminary to classify 
these objects as accreting magnetars. Note also that if at small X-ray luminosities convection is 
not developed in the shell, a quasi-Keplerian rotation law with n = 3/2 can be established. The 
equilibrium NS spin period in this case is P^ 1-3 ^ = Peg ^Rb/Ra ~ 10P e ^ -2 '. That means that 
NS periods in such XPSRs can easily reach a few 10000 s for the standard NS magnetic field. 

Other applications A possible implication of the settling accretion theory can be for non- 
stationary phenomena in XPSRs. The theory explains the observed temporary 'off-states in Vela 
X-l, GX 301-2, 4U 1907+09, when the plasma cooling near the magnetospheric equator occurs 
due to radiative processes [16]. The Compton cooling turns out to be ineffective due to X-ray 
pattern changing from fan-beam to the pencil-beam (as suggested by the observed X-ray pulse 
shape changes during the off-state in Vela X-l JT7[]). 

A dynamical instability of the shell on the time scale of the order of the free-fall time from the 
magnetosphere can appear due to increased Compton cooling and hence increased mass accretion 
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rate in the shell, leading to an X-ray outburst with duration lasting about the free-fall time scale of 
the entire shell (~ 1000 s). Such a transient behaviour is observed in supergiant fast X-ray tran- 



sients (SFXTs) (see e.g. []18[]), in which slow X-ray pulsations are found (e.g. IGRJ16418-4532, 
P* « 1212 s 119]). The observed flaring behavior can be the manifestation of a Rayleigh-Taylor 
instability from the magnetospheric radius occurring on the dynamical time scale ~ Ra / V GM. 

4. Conclusion 

At X-ray luminosities < 4 x 10 36 erg/s wind-fed X-ray pulsars can be at the stage of subsonic 
settling accretion. In this regime, accretion rate onto NS is determined by the ability of plasma to 
enter magnetosphere via Rayliegh-Taylor instability. The angular momentum can be transferred 
through the quasi-static shell via large-scale convective motions initiating turbulence cascade. The 
theory explains long-term spin-down in wind- fed accreting pulsars and properties of short-term 
torque-luminosity correlations. Long-period low-luminosity X-ray pulsars are most likely experi- 
encing settling accretion too. Spectral and timing measurements of slowly rotating X-ray pulsars 
can be used to further test this accretion regime. 
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